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ABSTRACT 

The current Swift sample of gamma-ray bursts (GRBs) with measured redshifts allows 
to test the assumption that GRBs trace the star formation in the Universe. Some 
authors have claimed that the rate of GRBs increases with cosmic redshift faster than 
the star formation rate, whose cause is not known yet. In this paper, I investigate the 
possibility for interpreting the observed discrepancy between the GRB rate history 
and the star formation rate history by the cosmic metallicity evolution, motivated 
by the observation that the cosmic metallicity evolves with redshift and GRBs prefer 
to occur in low metallicity galaxies. First, I derive a star formation history up to 
redshift z = 7.4 from an updated sample of star formation rate densities obtained by 
adding the new UV measurements of Bouwens et al. and the new UV and infrared 
measurements of Reddy et al. to the existing sample compiled by Hopkins & Beacom. 
Then, adopting a simple model for the relation between the GRB production and 
the cosmic metallicity history as proposed by Langer & Norman, I show that the 
observed redshift distribution of the Swift GRBs can be reproduced with a fairly good 
accuracy. Although the results are limited by the small size of the GRB sample and 
the poorly understood selection biases in detection and localization of GRBs and in 
redshift determination, they suggest that GRBs trace both the star formation and 
the metallicity evolution. If the star formation history can be accurately measured 
with other approaches, which is presumably achievable in the near future, it will be 
possible to determine the cosmic metallicity evolution with the study on the redshift 
distribution of GRBs. 
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1 INTRODUCTION 

Since the discovery of the afterglows of gamma-ray 
bursts (GRBs) and the determination of t heir redshifts 
l|Metzger et al.lll997l : Ivan Paradiis et al.lll997h , it has been 
firmly established tha t GRBs are at cosmological distances 
(for r e cent reviews see Piran et al.ll20oj : IZhang fc Meszaro3 
l2004l : iMeszarod bOod T" Observations on the hosts of 
GRBs have revealed that long-duration GRBs (hereafter 
GRBs) are associated with faint, blue and often ir- 
regular galaxies with high star formation rates (SFRs) 
llConselice et al.ll2005l: iFruchter et al.ll2006l: IXanvir fc LevanI 



Le Floc'h et al. 2006) . The disco ver y of the GRB-supernova 
conn ection (|Galama et al ] Il998l : Eil booel : IWooslev fc Heged 
|2006| . and references therein) supports the coUapsar model 
for long-duration GRBs (MacFadyen fc Woosley 1999; Mac- 
Fadyen, Woosley fc Heger 2001). 



Because of their 
be detected out to 
verse with minimal 
and d ust fciardi 
iBromm fc Loch 2002 



very high luminosity, GRBs can 
the edge of the visible Uni- 
extinction by intervening gases 



Loeb" '2000' 



Lamb fc RcichartI l2000l : 

Naoz fc Brombcrg 2007]) and 



12007 : Wainwright. Berger fc Penprasd 2007 . and references 
therein), confirming the early speculation that GRBs occur 
in star-formation regions and arise from the death of mas- 
sive stars (Paczyriski 1998; Wijers et al. 1998; see, however, 
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hence an ideal tool for probing the formation rate and 
the environments of stars at high redshift, the reion- 



tion jFynbo et al.' '2006a : Price et al 


'2006"; 'Savaglio' '2006|: 


Totani et al. 


2006; 


Bromm fc Loeb 


2007; Campana ct al.l 


20071: iGallerani et al.ll2007l: iProchaska et al.ll2007r). The ad- 



vantage of GRBs over quasars fo r probing the high reds hift 
Universe has been discussed by iBromm fc Loebl (|2007l ). It 
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found that flpRB/s fb. increases with redshift siRnificantly 
Da igne et ai]i20071: iLe fc Dermerll2007l: iKistler et al.ll2008l : 



Figure 1. The observed ratio of the GRB rate to the SFR, 
^GRB/SFRi ^ function of Q{z) and z, where Q{z) is defined as 
an increasing function of redshift z (equation 1171 and Fig. |6ll. Nor- 
malization is chosen so that ^Jgrb/SFR = 1 at Q = 0.5. The GRB 
rate history was obtained from a sample of b right Swift GRBs so 
that selection effects are minimized for 2^4 jKistler et al. lfioO^. 
see Section l6.ll of this paper). The SFR is the best fit to an up- 
dated sample of measured cosmic SFRs up to z ~ 7.4 (eqs.[T]and 
[2] Section [2]l. The horizontal dotted line denotes ^Jqrb/SFR = li 
an expected result of the assumption that GRBs trace the star 
formation unbiasedly. 



has also been proposed that GRBs can be use d as standard 
cand l es to constrain cosmoloKic a l parameters ijBloom et al.l 
I2OO3I : iFriedman fc BloomI l2005l : ISchaefeij 1200^ and refer- 
ences therein). However, this proposal has been seriously 
challenged by a recent study oflLil (|2007h . 

To study the relation between GRBs and the 
star formation, people often assume that the GRB 
rate is proportional to the SFR then compare the 
predicted distribution of the GRB redshift (or other 
parameters , e.g. the i n tensity) to the o bserved dis- 
tribution (iTotanil Il997l: iMao fc Mol Il998l: IWiiers et all 



Porciani fc Madaul I2OOII: iN ataraian et all l2005l: 



Jakobs son et al.l 20061 : Daigne et al. 2007; Le fc Dermerl 



20071 . for a review see Coward 2007) . If GRBs trace the star 



formation in the Universe unbiasedly, one would expect that 
the ratio of the GRB rate to the SFR (iicRB/SFR) does not 
vary with redshift. Then, an accurate measurement of the 
GRB rate history at high redshift would directly lead to the 
star formation history (SFH) in the early epoch which is 
otherwise hard to measure with the current technology be- 
cause of the uncertainty in dust obscuration for UV photons 
l|Hopkins fc Beacomlliooel ). 

Unfortunately, recent studies show that GRBs do not 
seem to trace the star formation unbiasedly (Fig. [TJ . Based 
on the current understanding on the SFH and the Swift 
sample of GRBs with measured redshifts, people have 



Yiik sel fc Kistleij l2007l : ICen fc Fanell2008l ). While observa- 
tions consistently show that the comoving rate density of 
star formation is nearly co nstant in the interval 1 < ^ < 4 
(IHopkins fc BeacomI |2006| ). the comoving rate density of 
GRBs appears evolving distinctly. 

Adopting a model-independent approach by selecting 
bright Swtft GRBs with the isot r opic-e quivalent luminosity 
Liso > lO^^erg s~ MKistler et al.1 (|2008l ) found that there are 
~ 4 times as many GRBs at redshift z ~ 4 than expected 
from star formation measurements. They claimed that some 
unknown mechanism is leading to an enhancement in the 
observed rate of high-redshif t GRB s. With a more sophisti- 
cated method, iDaigne et al.l l|2007l ) found that to reconcile 
the observed GRB redshift distribution with the measured 
SFH, the efficiency of GRB production by massive stars 
would be nearly six to seven times higher at 2 ~ 7 than 
at 2: ~ 2. Based on their results, Daigne et al. concluded 
that GRB properties or progenitors must evolve with cos- 
mic redshift. 

In this paper, I investigate the relation between the 
SFH and the GRB rate history, and explore the possibil- 
ity that the observed enhancement in the GRB rate at high 
redshift is caused by the cosmic metallicity evolution. Al- 
though the possibility of leading to an enhancement in the 
observed GR B rate evolution by the cosmic metallicity was 
mentioned bv lKistler et al.l l|2008h . they did not give a quan- 
titative analysis or a detailed discussion. Instead, Kistler et 
al. discussed more thoroughly on other possible causes, in- 
cluding evolution in the fraction of binary systems which 
had been proposed as a channel for producing GRBs, an 
evolving initial mass function (IMF) of stars, and evolution 
in the galaxy luminosity function (LF). 

There is growing evidence that metallicities play an im- 
portant role in the production of GRBs. Observations on 
the hosts of GRBs revealed that GRBs prefer to occur in 
galaxies with low metallicities (Fynbo et al. 2003, 2006a; 
Prochaska et al. 2004; Soderberg et al. 2004; Gorosabel et al. 
2005; Berger et al. 2006; Savaglio 2006; Stanek et al. 2006; 
Wolf & Podsiadlowski 2007; Modjaz et al. 2008; Savaglio, 
Glazebrook fc Le Borgne 2008). Based on the observa- 
tional and theoretical evidence that the mass-loss rate of 
Wolf-Rayet stars depends on the metallicity ( Crowther et al.l 
l2002l : IVink fc de Koted l2005l : ICrowtheil 12007^), theoretical 
studies on the coUapsar model of GRBs arising from sin- 
gle massive stars suggested that GRBs can only be pro- 
duced by stars with metallicity Z < O.IZq since otherwise 
strong stellar winds will cause stars to lose too much mass 
and angular momentum to form a disk around a black hole 
of several solar masses which is essential for the produc- 
tion of GRBs (Hirschi, Meynet & Maeder 2005; Yoon & 
Langer 2005; Woosley & Heger 2006; Yoon, Langer fc Nor- 
man 2006). 

It is well-known that the cosmic metallicity evolves 
strongly with redshift, and galaxies at higher redshift tend to 
have lower metallicities (Pettini et al. 1999; Prochaska et al. 
2003; Rao et al. 2003; Kobulnicky fc Kewley 2004; Kewley & 
Kobulnicky 2005, 2007; Kulkarni et al. 2005, 2007; Savaglio 
et al. 2005, 2008; Wolfe, Gawiser fc Prochaska 20 05; Savaglio 
2006; Peroux et al. 2007). INataraian et all (|2005l ) considered 
a model where GRBs trace the average metallicity in the 
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Universe rather than the SFR, and the GRB rate decreases 
with increasing metallicity. However, their model led to a 
GRB redshift distribution that is nearly indistinguishable 
from the distribution predicted by the SFR (jjakobsson et al] 
I2OO6I ). 

Recently, iLanger fc NormanI (|2006l ) considered the ef- 
fect of the cosmic metallicity evolution on the integrated 
production rate of GRBs in the framework of the coUap- 
sar model, assuming that the GRB rate is jointly deter- 
mined by the SFR and the metallicity evolution. Adopt- 
ing a simple model for the metallicity evolution and a best 
fitted SFR, and assuming that a GRB is produced if a pro- 
genitor star is massive enough and has a metallicity below 
a threshold {Z < O.IZq), they showed that the observed 
global ratio of the GRB rate to the core-collaps e supe rnova 
rate (~ 0.001) can be reproduced. iNuza et al.l (|2003) and 
ICen fc Fand (I2OO8I ) investigated the host galaxies of GRBs 
in a cosmological hierarchical scenario with numerical sim- 
ulations. They found that the observed properties of GRB 
hosts are reproduced if GRBs are required to be generated 

by low metallicity stars. 

I will incorporate the model of iLanger fc NormanI 

l|2006ll into a probability distribution function of the lumi- 
nosity and redshift of GRBs to study the rate history of 
GRBs. For this purpose, I will present an updated SFH ob- 
tained by adding th e new measuremen ts of SFR densities at 
z ~ 2.3 and 3. 05 bvlReddv e t al. (2008) and at z ~ 3.8, 5.0, 
5.9 and 7.4 bvlBouwens et al (2007 . 2008) to the data com- 
piled bv lHopkins fc Beacom ( 20061 ). With the updated data 
I will derive an analytic formula for the SFH and show that 
the observed distribution of Swift GRBs can be successfully 
reproduced when the evolution of the cosmic metallicity is 
properly taken into account. Then I will argue that, after 
a significantly expanded and well-defined sample of GRBs 
with measured redshifts and luminosities is available in fu- 
ture and the SFH at high redshift is accurately determined 
with other approaches, GRBs will be a powerful tool for 
probing the cosmic metallicity evolution. 

The paper is organized as follows. In Section[21 1 present 
the updated measurements on the SFR up to z ~ 7.4 and 
derive an analytic formula for the cosmic SFH by fitting 
the updated data. In Section [31 I summarize the current 
measurements on the evolution of the metallicity in galaxies 
with cosmic time and argue that the data indicate a consis- 
tent picture for the cosmic metallicity evolution when the 
redshift-dependent relation between the metallicity and the 
galajcy stellar mass is considered. In Section|4l I describe the 
Swift GRB sample that is used for the current work, show 
the luminosity distribution of the GRBs, and discuss the se- 
lection biases involved in the detection and localization of 
Swift GRBs and the measurement of their redshifts. In Sec- 
tion O the model that I adopt for calculating the GRB rate 
history is outlined, which includes assumptions about the 
probability distribution function, the SFR, the evolution of 
cosmic metallicity, and the form of the GRB LF. In Section 
[6] I present the results calculated with the model, and fit the 
observed distribution of the luminosity and redshift for the 
whole GRB sample and a bright GRB subsample. In Section 
[T] I draw conclusions and discuss some implications of this 
work. 

Throughout the paper, I assume a flat universe with 
fim = 0.3, = 0.7, and Ho = 70 km s"^ Mpc"\ 



2 THE UPDATED STAR FORMATION RATE 
HISTORY 

With modern UV and far-infrared (FIR) observations the 
SFH has been well established fo r redshift z < 4, with espe - 
cially tight constraints for z < l. lHopkins fc BeacomI (|2006l ) 
compiled and critically analysed the data that were available 
then. They found that, despite large data scatter, the SFR in 
the redshift range of 1 < z < 4 is approximate ly a constant, 
which agreed with the previous claim (e.g., ISteidel et al.l 
Il999l ). Beyond z ~ 4, although the data were highly incom- 
plete, a meaningful constraint on the SFH was drawn: the 
SFR declines with z > 4. From the 120 data points collected 
from UV, FIR, radi o, Hq and the Hubble Ult ra Deep Field 
(HUDF) estimates, iHopkins fc BeacomI (|2006l ) critically se- 
lected 56 'good' data points and fltted the SFH with simple 
analytical formulas and derived conservative uncertainties Q 
For the criteria applied in the selection of 'good' data 
and a complete lis t of re f erence s for the data, please re- 
fer to Hopkins & Beacornl (|2006|) . Here I update the sam- 
ple of lHopkins fc Beacom (2006) by including the new mea- 
surern ents of iBouwens et al.. (2007. . ,2008. ) and iReddv et al.l 
l|2008l ). 

Using the HUDF and the GOODS fields. iBouwens et all 
(2007, 200S) found large samples of star- forming galaxies at 
z ~ 4, 5, 6, and 7-10. With those data, the rest-frame UV 
LFs were determined with high accuracies at z ~ 4-6. It has 
been found that the faint-end slope a and the normalization 
factor (p* show very little evolution with cosmic time, but 
the characteristic absolute magnitude M/jy brightens con- 
sider ably from 2 ~ 6 to ~ 4 (by ~ 0.7 mag) (jBouwens et al.l 
|2003). With all available deep optical and near-IR (NIR) 
data over the two GOODS fields. iBouwens et all (|2008l ) have 
also derived a rest-frame UV LF at z ~ 7.4, and obtained a 
constraint on the UV LF at z ~ 10. The SFR densities at 
z ~ 4-10 were derived from those UV LFs, confirming that 
the cosmic SFR density decreases quickly with increasing 
redshift beyond z ~ 4. 

Using a large sample of rest-frame UV-selected and 
spectroscopically observed galaxies in the redshift interval 
1.9 ^ z ^ 3.4 combined with ground-based spectroscopic 
Hq and Spitzer MIPS 24 fim data, which includes over 
2000 spectroscopic redshifts and ~ 15000 photometric can- 
didates in 29 independ ent fields covering a total area of al- 
most a square degree, iReddv et al] (|2008l ) derived robust 
measurements of the rest-frame UV, Ha, and IR LFs at 
1.9 < z ^ 3.4. The results indicate that the UV LF under- 
goes little evolution between 2 ~ 4 and z; ~ 2. The SFR 
density at z = 2.3 ± 0.4 and z = 3.05 ± 0.35 was derived 
from the UV luminosity de nsity and the IR l uminosity den- 
sity respectiv ely, using the KennicuttI (|l998l ) relations and 



sity respectiv ely, using tne [ ivenmcuttl l|iijUel l relations 
assuming the lSalpete"T(|l955l ) IMF from 0.1 to 100 Mq . 

To include t he new data poin ts of iBouwens et al.l 
(2007, 20 0^) and iReddv et al.i (|2008i ) into the sample of 
iHopkins fc BeacomI (|2006l ). two tvpes of corrections must be 
considered. First, UV lights are strongly obscured by dust 
so a dust-obscuration correction factor Ci must be applied 
to the SFR density derived directly from the UV luminosity 



^ The number 5 8 of 'good' data points printed in 
iHopkins &: BeacomI | j2006l) was a typo (A. M. Hopkins, pri- 
vate communications). 
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Figure 2. Dust corrections for UV continuum emissions at dif- 
ferent redshift intervals. The hatched r egion is the du s t cor - 
rection obtained from the FIR SFH of iLe Floc'h et alj ||2005|') 

1z < 1) and the FIR measurements of iPerez^GonzHezer'al] 
20051 ) (1 < 2 < 3) (see text and H opkins &: Beacom 2006 fo r 
details), relative to the UV SFH of ISchiminovich et al.l l l2005h . 
The two pentagons are estimated from the IR and UV SFRs of 
iReddv et al]||2008h . Ci 5.17 at z ~ 2.3 and Ci 4.55 at z ~ 3. 
T he triangle is the d ust correction Ci ~ 3.80 at 2 --^ 4 derived 

by 
by 



' OuchietaL [f2004). The circle is the dust correction estimated 
Bouwens et al.. (^006, 2002), Ci ^ 1.51 at 2 ~ 6. The dashed 
lines are the linear interpolation between 2 = 4 and 2 = 6, and 
extension beyond 2 = 6. Note, the dust corrections beyond 2 ~ 4 
are highly uncertain. 



density. For z ^ 1, the FIR SFR de nsit y was well measured 
with S pitzer (|Le Floc'h et aljlioosh . In iHopkins fc Beacoml 
(|2006l ). the UV data at z < 1 were 'obscuration cor- 



rected' by adding the FIR SFR density from lLe Floc'h et al.l 
l|2005l ) to each point. For obscuration corrections to the 
UV data between z ~ 1 and z ~ 3, iHopk ins fc Beacoml 
(l2006l) made use of the fact that the FIR measurements of 
iPerez-Gonzalez et al are quite fiat i n this domain as 

well a s being highly consistent with those of lLe Floc'h et all 
l|2005l ) at z < 1, an d add the constant SF R density corre- 
sponding to that of iLe Floc'h et al.l ()2005l ) at z = 1. This 
looks a reasonable treatment for the obscuration correction 
at 1 < z < 3 since later measurements of the IR luminosity 
density (jCaputi et al.1 l2007l: iReddv et al.1 12OO8I) ag ree with 
a the trend found by Perez - Gonzalez et al.' !2003) in this 
redshift range (fig. 27 of iReddy ct al. 2008). In addition, 
the dust cor r ection s derived from the IR and UV SFRs of 
iReddv et al.1 (|2008l ) are Ci « 5.17 at z ~ 2.3 and Ci ^ 4.55 
at z ~ 3, in good agreement with the values of 4.5-5.0 at 
z = 1.0-3.5 o btained with other approaches (ISteidel et al.l 
ll999l : lNandra ct al. 2002; Ro ddy fc Steidel 20^^ 

Beyond redshift z ~ 3, dust obscuration is quite un- 
certain. Investigating the photometric properties of Ly- 
man br eak galaxies (LGB s) detected in Subaru Deep Field 
(SDF). IOuchi et al.l (|2004h estimated the dust extinction of 
z ~ 4 LGBs and obtai ned E{B - V) = 0.15 ± 0.03. By 
the lCalzetti et al] l|2000l ) extinction law, this redenning cor- 
responds to a dust correction factor Ci « 3.80. The dust 
extinction at z ~ 6 is still unknown. Observing the fact 



that the UV continui um slope /3 at z ~ 6 is bluer than 
that observed at z ~ 3 JStanway. McMahon fc Bunkerl 



Yan et al.|[2005l : lBouwens et al.li2006h . lBouwens et al 



2005; 



20061 , 



2007^ argued that the dust extinction at z ~ 6 is lower than 
that at z ~ 3. Adopting /3 = —2 (|Bouwens et al.ii2006.) and 
the e x tinction-/3 relati o n Aiggo = 4 .43-1-1.99/3 ( Meurer et al.l 
ll999l) . lBouwens et al.l (|2007l . I2OO8I ) assumed that Ci « 1.51 
(0.18 dex) at z ~ 6 in their SFR derivation. However, com- 
parison of the SFH derived from the Local Group and other 
nearby galaxies with the cosmic SFH indicates a factor of 
~ 10 extinction correction to high-redshift and UV-based 



SFR measures (|Drozdovskv et al.|[200i ). 

I wiU adopt Ci « 3.80 at z ^ 4, and Ci « 1. 51 at 
z > 6. At 4 < z < 6, I wiU follow iBouwens et"all (|2007| . 
l2008h and estimate the dust correction by linear interpola- 
tion between the z ~ 4 and z ~ 6 values. The treatment 
for dust corrections outlined above is sketchedly presented 
in Fig. H 

The next correction that should be considered is the 
factor that converts a 'partial' SFR calculated by integra- 
tion down to a limit luminosity Lmin > to a 'total' SFR 
calculated by integration down to L = 0, which is simply 
given by C2 = r(2 -I- a)/r(2 + a, Lmin/L*), where a is the 
faint-end slope of the Schechter LF, and L* is the charac- 
teristic luminosity. Here r(a:) is the gamma function, and 
r{a,x) is the incomplete gamma function. 

The 'partial' SFR calculated down to the survey mag- 
nitude limit is the SFR often reported in the literature. The 
'total' SFR is strongly denpendent on uncertainties in the 
faint-end slope, but it allows direct comparison with differ- 
ent SFR measurements (Schiminovich ct al. 2 0051). In th e 
compilation of the SFR data in lHopkins fc Beacoml (|2006l ). 
it has been attempted to integrate to L = in all cases 
where possible, except for several data points which will 
be explained and remedied below (A. M. Hopkins, private 
communications). To add the new data to the sample of 
iHopk ins fc Beacoml l|2006l) . the integration correction must 
be made whenever it is relevant. 

In the calculation of the integration c orrection factor 
C2, I make use of the UV LFs derived by H ouwens et al.l 
(|2007l . |2008| ) in the redshift interval 4 < z < 10. Since the 
fain-end slope a of the LFs shows very little evolution with 
z, a mean value of —1.71 is used in all cal culat ions. 

The new data of lBouwens et al.1 (|2008h and lReddv et al.l 
(200^, with the above corrections being applied, are 
added to the 'good' data compiled by iHopkins fc Beacoml 
(20 03) and shown in F ig. [3] Several data points of 
Hopkins fc Beacom ('2006') have been modified or removed. 
Since Bouwens ot al. (2006)'s SFR at z sa 6 has been up- 
dated to the new measureme nt by Bouwens et ahl (|2007| . 
i2008i ). it has been removed from lHopkins fc Beacoml (l2006l )'s 
sa mple (then the sample c ontains 55 data points). The SFRs 
of'Giavalisco et al.' (|2004h at z « 3.78 and 4.92 and that of 
[Bunker ct al. (2004) at z « 5.9, all in iHopkins fc Boaco^ 
(2006)'s sample, were only integrated to a finite Lmin (= 
0.21/^^3 and O.lL'^a, respectively; A. M. Hopkins, private 
communications). They are updated by including the inte- 
gration correction and with new dust correctio ns according 
to the procedure described above. The SFR of lOuchi et al.l 
(|2004l ) at z ~ 4.7 is updated with a new dust correction 
factor obtained by linear interpolation between z = 4 and 
z = 6. 
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Figure 3 . Updated star f o rmati on history. The sample contains 51 data points from iHopkins fc BeacomI 1I2OO6I} (blue) , two UV data 
points of iGiavaUsco et al.l ||2004|) with updat ed corrections (red squares), one UV data point of iBunkereT'an l l2004i with updated 
corrections ( r ed sta r), one UV data point of lOuchi et al.l l|2004h with updated corrections (red triangle), two UV+IR data points of 
iReddv et al.l 1 I2OO8I) (red pentagons), four UV data points and an upper limit of iBouwens et all {2008) (red circles and the downward 
arrow) (see text for details). The total number of data points (not including the upper limit) is therefore 61. The solid black line is the 
best fit wi th a piecew i se pow er-law (eqs. [T]and[2]l, with Xv — 1-18 (55 degrees of freedom). The dashed black line is the best fit with a 
formula of lCole et al.l JgOOy) (eq. [Sjl, with = 2.44 (57 degrees of freedom). In the calculation of chi-squares the errors in 2 (bin sizes) 
are not taken into account. 



Table 1. The star formation rate density derived fro m the new UV observational da ta llBouwens et al.ll200^ ) and from the old UV 
data with new corrections (Bunker et al. 2004; Giava lisco et al .^ 2004; Ouchi et al. 2 0041 ) 



Reference Estimator Redshift logSgpRun" logSsFR* logCi'^ logC2'' 



Giavahsco et al. (2004) ~ 1500 A 3.780 ± 0.340 -1.639 ±0.060 -0.562 ± 0.060 0.580 0.496 

4.920 ±0.330 -1.855 ±0.140 -0.902 ±0.140 0.379 0.575 

Bunker et al. (2004) i-Dropouts 5.900 ± 0.500 -2.301 ± 0.079 -1.627 ±0.079 0.180 0.494 

Ouchi et al. (2004) - 1500 A 4.700 ± 0.500 -1.521 ±0.079 -1.097 ±0.079 0.424 0.000 

Bouwens et al. (2008 ) B-Dropouts 3.800 ± 0.300 -1.720 ±0.050 -0.644 ± 0.050 0.580 0.496 

y-Dropouts 5.000 ±0.500 -2.050 ± 0.060 -1.112 ± 0.060 0.363 0.575 

i-Dropouts 5.900 ±0.300 -2.180 ±0.080 -1.311 ± 0.080 0.180 0.689 

z-Dropouts 7.400 ±0.700 -2.580 ± 0.210 -1.548 ± 0.210 0.180 0.852 

J-Dropouts 10.00 ±1.200 < -2.760 < -1.638 0.180 0.942t 



"Uncorrected SFR density SgpH in units of Mq yr ^ Mpc ^. 

''Dust and integration corrected SFR density, Sgpj)^ = C iC2SsFR.uni iri un its of Mq yr~^ Mpc~'^. 

''Dus t correction factor Ci. At 2 ~ 4, I adopt Ci = 3.80 l lOuchi et al.ll2004h . At 2 > 6, I adopt Ci = lO'''^* = 1.51 llBouwens et al.l 
[2003). For 4 < 2 < 6, the value of C\ is obtained by linear interpolation between 2 = 4 and 2 = 6 (Fig. [2]l. 

''Integral correction factor C2, converting a 'partial' SFR defi ned by integration dow n to luminosity L = Lniin to a 'total' SFR 
defined by integration do wn to L = 0. The UV LF s pu blished in [ Bouwens et al.l l|2008h are adopted, with an a verage faint-end slop e 
a = -1.71. For the data of lGiavalisco et al.l ||2004|) a nd lSouwens e t al. (200i), Lmin = 0.2L*^3. For the data of lBunker et al.l (i2004l) 
Lmin = 0.1-L*_3. For the data of lOuchTet al.l ||20M ). the SsFR,un has already been integrated down to L = bv lHopkjns'irBeacoml 
1I2OO6I ) so C2 = 1. 

^C*2: For the upper limit at 2 = 10 of iBouwens et aL I | |2008|) . I take the characteristic magnitude = —19.6 in the calculation 

of C2. 
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Table 2. lReddv et alj l|2008h 's star formation rate density derived 
from UV and IR measurements (in Mq yr~^ Mpc~'^) 



that at z > 4 the data suffer large uncertainties due to the 
lack of knowledge in dust obscuration as explained above. 



Redshift 


logSsFR.Uv" 


logSsFR.IH,'' 


logSsFR^ 




2.30 ±0.40 


-1 2g4+°-°'"5 


-0 674+° °'^* 


-0.58llHi^ 


3 THE COSMIC METALLICITY 


3.05 ±0.35 


-1 4Q4+0060 


-0 944+0-<"'5 




Observations have shown that the metallicity 



"Uncorrected SFR density derived from the UV luminosity den- 
sity, integrated down to L = 0.04L*_g. 

'SFR density derived from the IR luminosity density between 8 
and 1000 fim. 

'^The total SFR density obtained by the sum of SsFR UV 
SsFR,lR. 

In Fig. |3l all new data (6 points, not counting the 
upper limit) and updated old data (4 points) are shown 
wi th red symbols, and the unchanged data in the sample 
of lHopk ms fc BeacomI (|2006i ) are shown with blue symbols. 
The new data are broadly consistent with previous mea- 
surements, and extend the sample to higher redshifts (up to 
z ~ 7.4). The new data and the upda t ed old data are listed 
in Tables [Hand [2] Since iReddv et al.l (|2008h 's SFRs are de- 
termined by combination of UV and IR measurements, no 
extra dust correction is needed. The UV SFR is obtained by 
integration down to 0.04L*=3 which is already small enough. 
Since the UV SFR contributes to the total UV-I-IR SFR only 
by a small fraction (~ 20%), a further i ntegration cor r ection 
is not necessary. Hence the SFR data of lReddv et al.l (|2008l ) 
are listed in a separated table. 

In Fig. Eland Tables [1] and [2] the Salpeter IMF has been 
assumed. Conversion of SFH estimates to an alternative IMF 
ass umption correspond s to a s imple scale factor, as described 
bv lHopkins fc BeacomI l|2006l ). 

For the latter application, the SFH presented in Fig. 
5] is fitted by a piecew ise power-law, as was did in 



Hopkins fc BeacomI (|2006l ). The results are 



log Esfr(«) = a + 61og(l -I- z) 




(-1.70,3.30) , 
(-0.727,0.0549) 
(2.35, -4.46) , 



z < 0.993 
0.993 <z < 3.80 
z > 3.80 



(1) 



(2) 



and EsFR is in units of Mq yr~^ Mpc~"^. The reduced chi- 
square of the fit is x? ~ 1-18 (55 degrees of freedom; errors 
in z not counted). The value of Xr indicates that the fit is 



The SFH is also fitted with a formula of ICole et al.l 



a + bz 
1 + (z/cY 



(3) 



The results are {a,b,c,d) = (0.0157,0.118,3.23,4.66), with 
2.44 (57 degrees of freedom). 



Compared to the fit results of iHopkins fc BeacomI 

l|2006h . the updated sample leads to a fiatter SFH at 1 < 
2 < 4, and a less rapidly decreasing SFH at z > 4. A ma- 
jor reason for the large differenc e in the fi tted high-z slope 
(^ _4 versus ~ —8) is that in iHopkins^ BeacomI (|2006| ) 
the SFR density of Bunker et al. at z = 5.9 has been under- 
estimated by a factor ~ 3 since it has been integrated only 
down to Lmin = 0.11/^^3. However, it should be stressed 



in galaxies decrease with increasing redshift. This is in agree- 
ment with the hierarchical scenario of structure formation, 
in which metals are assumed to be produced by the forma- 
tion of stars, injected into the interstellar medium (ISM) by 
supernova explosions, and some of the metals are expelled 
from the host galaxy by supernova-driven galactic winds. 
However, measurements of the cosmic metallicity evolution 
with different approaches do not appear to converge, and 
numerical simulations of galaxy formation are still not suf- 
ficient for fully reproducing the observations on the cosmic 
metallicity evolution (Nagamine, Springel fc Hernquist 2004; 
Kobayashi, Springel fc White 2007). 

Metallicities in damped Lya absorbers (DLAs) in the 
spec tra of background QSOs have been exten si vely stud- 
ied dPettini et all 1 19991: IPrqchaska et a .1 l2003l : iRao et all 



I2OO3I : iKulkarni et all l2005l : IWolfe et all \200^ . The data 



consistently suggest that the metallicity in QSO-DLAs 
evolves with cosmic time according t o ZIZp, oc 10 . with 
7 K. 0.2-0.4 JProchaska et al.ll2003l: iNagamine et al.l 12004 : 



Kulkarni et al.1 |2005| . l2007l : ISavaghd I2OO6I : iPeroux et all 
2007h . In addition, there is evidence indicating that the 



metallicity in sub-DLAs (absorbers with H I column density 
A^'hi < 10^"'^ cm~^) evolves with redshift more strongly than 
that in normal DLAs (with A^hi > 10'^°'^ cm~^), and a re 
more metal-rich (|Kulkarni et al.ll2007l : IPeroux et al.ll2007l '). 

DLAs have also been observed in the afterglows of 
GRBs, indicating metallicities and column densit ies that 
are generally higher than in QSO-DLAs (jVreeswiik et al-l 
2004l:lBerger et al.ll2006l : |Prochter et al.ll2006l:ISavagholl2006l; 



Sudilovskv et al.l l2007l : IWatson et aP l2007l: iFynbo et all 
2OO8I I. A system study on GRB-DLAs l|Savaglioll2006l ) in- 



dicates that the metallicity in GRB-DLAs evolves with red- 
shift with a slower rate (7 ~ 0.18) than that in QSO-DLAs, 
suggesting that GRB-DLAs and QSO-DLAs belong to two 
different classes of absorbers although they both trace the 
ISM in high-redshift galaxies. 

By measuring nebular oxygen abundances for 204 
emission-line galaxies with redshifts 0.3 < z < 1.0 in the 
Great O bservatories Origins Deep Su rvey-North (GOODS- 
N) field. iKobulnickv fc KewlevI l|2004h found that the metal- 
licity in galaxies with —18.5 < Mb < —21.5 evolves with 
7 ~ 0.14 from z = to 1. Compilation of oxygen abun- 
dances in star-forming galaxies with AIb < —20.5 and 
< z < 3.5 indicates that the metallicity in star-forming 
galaxies evol ves with 7 » .15 in the redshift range of 
< z < 3 tKewlev fc Kobulnickv 2005. 2007). 

Different types of metallicity measurements suffer dif- 
ferent selection biases so the fact that the values of 7 do 
not converge is probably not surprising. Detection of ab- 
sorption lines in the neutral ISM of galaxies acrossing QSO 
sight-lines gives information of one line of sight in the galaxy 
so the result suffers large statistical fiuctuations. Detecting 
emission lines in the integrated galaxy spectra is observa- 
tionally challenging at high redshift since prominent lines get 
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very weak and redshifted to the NIR range (|Savaglio et alj 
120051 ). The different evolution observed for QSO-DLAs and 
GRB-DLAs may indicate that GRB-DLAs are associated 
with more massive g alaxies, on average as massive as the 
LMC llSavaglid 12009 ). Reconciling the metallicity distribu- 
tion of GRB-DLAs, QSO-DLAs, and Lyman-break galaxies 
(star- formin g galaxies) at high redshift has recently been 
discussed by jPynbo etal] (|200^ ). 

It is well-known that the stellar mass of star-forming 
galaxies is correlated with the metallicity: galaxies with 
larger stellar masses tend to have higher metal li cities (the 
mass-mctallicity relation, 'Sava glio et al.l l2005l : lErb et al] 
2006; Kowlcy & Ellison 2008, and references therein). 
Savaglio et al.l l|2005l ) derived a simple relation that gives 



the metallicity in the gas of a galaxy at a given ste llar mass 
and H ubble time (their eq. 11). With this relation. ISavagliol 
l|2006 t) argued that the apparent difference in the values of 
7 for QSO-DLAs and GRB-DLAs may be explained by the 
mass-metallicity relation with the assumption that different 
measurements probe galaxies with different stellar masses. 

Using the scaling 12 + log(0/H) = \og{Z/ZQ) + 8.69 
l|Savaglidl20"0^ . Savaglio's relation can be written as 



\og{Z/ZQ) = -16.28 2.53 log M- 0.0965 log^M 
-h5.171ogte - 0.394 log^te 
-0.403 log log A/ , 



(4) 



where M is the stellar mass in units of Mq, and tn is the 
Hubble time in units of Gyr. 

The relation was used to fit the evolution of the 
metallicity in QSO-DLAs and GRB-DLAs. It was pre- 
dicted that the stella r mass of QSO-DLAs is lO^'^^^^ ^'^Afo 
ISavaglio et al.ll2005l ). and the stellar mass of GRB hosts is 
preferably in the range of 10*'^-10''-*AfQ l| Savaglio! [200^ . 
The latter is in agreement with the average stellar mass 
of lO^'^Af0 with a 1-a dispersion of 0.8 dex found for 46 
GRB hosts with opt ical-NIR photometry and spectroscopy 
jSavaglio et al.ll2008D . 

With a galaxy stellar-mass function ^{M), two aver- 
age metallicities can be defined. The first one is defined by 
averaging over the number of galaxies 



Z 

Z^ 



J^^^^ Z<f{M)dM 



(5) 



For a Schechter function of galaxy stellar-masses (Panter, 
Heavens & Jimenez 2004), the low- mass-end slope a is usu- 
ally < —1. Then, the integral of <E>(Af) over M does not con- 
verge as the lower limit of the integral goes to zero. Hence, 
a non-zero minimum stellar mass Mmin is necessary in the 
definition of {Z/Zq)i. 

The other average metallicity is defined by averaging 
over the stellar mass (or the mass of H I gases) 



Z 



ZM${M)dM 
Zq M^(M)dM 



(6) 



i.e., the total metal mass divided by the total stellar (or Hi) 
mass. In the definition of {Z/Zq)2, the lower limit of mass 
in the integral is set to zero since the integral converges. 

Adopting a Schechte r function of galaxy stellar-masses 
from IPanter et all (|2004l ) with a = —1.16 and the critical 
mass Af* = 7.64 x 10^°h~'^MQ (with h = 0.7, h is the 



SI 

\ 




Figure 4. The cosmic metallicity as a function of the redshift. 
The solid line is the fit to the metallicity in star-forming galax- 
ies by Z/Z (7) oc with 7 = 0.15 l lKewIev fc Kobulnickvl 
I2OO5I.I2OO7I) . The two dotted lines are the fits to QSO-DLAs with 
7 = 0.352 (low er line) and GRB-DLAs with 7 = 0.18 (upper line; 
ISavaelidlioO^) . The two dashed lines are the average metallic i- 
ties calculated with the empirical model of lSavaglio et all ll2005t) . 
{Z/Zq)i (lower line, eq. [5j and (Z/Zq)2 (upper line, eq. O. 



Hubble constant in units of 100 km s ^ Mpc ^), I calcu- 
lated (Z/Zq)i (with Afmin = lO'^-^A^o) and {Z/Zq)2 for the 
metallicity given by equation (|4|. The results, with compar- 
ison to measurements, are shown in Pig. 3] 

The {Z/Zq)i defined by equation (O is suitable for 
comparison with observations since the observed metallic- 
ity evolution is usually derived from the fit to the distribu- 
tion of metallicity abundances in individual galaxies. Prom 
Fig- 131 {Z/Zq)i (the lower dashed line) is between the mea- 
sured metallicity in star- forming galaxies (with 7 = 0.15, 
solid line) and that in GRB-DLAs (with 7 = 0.18, upper 
dotted line). 

The {Z/Zq)2 defined by equation ((6)l (or its variant, 
e.g., ratio of the total mass of metals to the total mass of H I 
gases in all ga laxies) is sometimes ad opted in numerical sim- 
ulations (e.g.. lNagamine et al.ll2004h . Prom Pig.H {Z/Zq)2 
(the upper dashed line) gives an estimate of the metallic- 
ity higher than that derived from the fit to measurements, 
and evolves with z with a slower rate. This is caused by 
the fact that, by its definition, the contribution to {Z/Zq)2 
comes dominantly from galaxies with stellar masses around 
Af* ~ 10^^ Mq while faster evolution and lower metallici- 
ties are attributed t o galaxies with smaller stellar masses 
ISavaglio et al.ll2005l ). 

Por a < —1, the value of {Z/Zq)i and its evolution are 
sensitive to Afmin, since the dominant contribution to the 
integral comes from low stellar-mass and hence low metal- 
licity galaxies because of their large numbers. The choice of 
seems reasonable according to the mea- 
surement on_thB_distriMition of stellar masses for GRB host 
galaxies ijSavaglio et al. 

Em,)- 
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Although SavagUo's relation leads to reasonable results 
in fitting the metallicity distribution of QSO-DLAs and 
GRB-DLAs by varying the stellar mass, its validity and cal- 
ibration need to be tested and improved with future mea- 
surements on metallicity and stellar masses of galaxies. For 
the purpose of fitting the GRB rate history in this paper, it 
is enough to adopt an empirical evolution equation 

Z/Zq oc 10"^" (7) 

with a constant 7 as adopted bv lLanger fc NormanI l|2006f ). 
Given the limit number of GRBs with measured redshifts 
and spectra and many unknown observational biases, it is 
impossible to constrain the parameters in the metallicity 
evolution with GRBs at the present time. 

Since long-du ration GRBs occurred in star-forming 
galaxies, following iLanger fc NormanI 1I2OO6I) I will adopt 
the v alue 7 — 0.15 derived by Kewlev fc Kobulnickvl l|2005l, 
l2007l 'l. 



4 THE SWIFT GRB SAMPLE 

4.1 The Luminosity Distribution of Swift GRBs 

Since the launch of Smft in late 2004, 310 GRBs have 
been detected by 31 March 2008, of which 181 bursts have 
optical detections (by the UV/Optical Telescope on Swift 
and/or ground telescopes) and 102 bursts have measured 
redshiftsl30 

JButlor ct al.! (|2007l ) have compiled a catalog of 218 Swift 
GRBs and calculated their durations and spectral parame- 
ters between and including GRBs 041220 and 070509, in- 
clu ding 77 events with measured redshifts. From that cata- 
log, iKistler et al.l (2OO8) selected 63 bursts with long dura- 
tions (T90 > 2 s, T90 is the observed time duration to contain 
90% of the total counts, with 5% in front and the other 5% 
behind.) and reliable redshift measurements to investigate 
the redshift distribution of Swift GRBs. The following four 
long GRBs are contained in the catalog of Butler et al. but 
not included in the sample of Kistler et al., because of their 
unrehable redshifts: 060116, 060202, 060428B and 061004 
(M. D. Kistler, p r ivate communications) . 

iLiang et al.f (|2007l l have derived the luminosity for 
45 long Sw i ft GR B s using the method developed by 
IZhang et all (|2007bl ). I Cabrera et"al] l|2007h have published 
the spectral and energetic results for a smaller sample of 
SwiftGK^s, which contains 29 long bursts with calculated 
isotropic-equivalent enerCT^ 

I use th e samp le of iKistler et al.l (|2008t ) selected from 
iButler et all l|2007l) — with addition of GRB 050223 —for 
my analysis, since it is the largest sample to date. The 
64 GRBs in the sample have determined durations, red- 
shifts, peak spectral energy, and the isotropic-equivalent 
energy in the 1-10'* keV band in the rest frame of 
GRBs (|Butler et al.ll2007l : iKocevski fc Butlej|2007h . For the 
isotropic-e quivalent energy -Ei a o, I ta ke the values from the 
table 1 of iKocevski fc Butler! (|2007l ) since they have more 
si gnificant digits than the corresponding values in the table 2 
of iButler etakl l|2007l ). However, GRB 050223 (not included 



^ http:/ /swift. gsfc.nasa.gov/docs/swift/archive/grb_tablc/ 
^ http://www.mpe.mpg.de/ jcg/grbgcn.html 




48 50 

l°g L,s<. (erg s-') 



54 



Figure 5. Distribution of the isotropic-equivalent luminosity (de- 
fined by eq. [Sj for 64 long-duration Swift GRBs (the solid line 
histogram, with the number of GRBs in each bin indicated by a 
dark point with Poisson error bars) . The solid curve is a Gaussian 
fit to the distribution of logLiso, with x? = 1.24. The dashed 
curve is a least chi-squares fit by the model in Section [S] with 



1.18 (see Sections [5] and [6] for details). The 45 bursts to 



the right of the vertical dashed line at E\ao = 0.8 X 10 erg s~ 
form a subsample of 'bright GRBs', which are not aff'ected by the 
luminosity threshold when z < A (Section 16. 111 . 



in Kis tler et al. 2008) was not listed in IKocevski fc Butlerl 
(l2007l). so I take the value of Eiso for this event from 
iButler et all (|2007l). 

Following IKistler et al.l ([2OO8), I calculate the isotropic- 
equivalent luminosity of a GRB by 



Eiso 
Tgo 



(1 + ^) 



(8) 



The distribution of Liso for the 64 GRBs in the sample is 
shown in Fig. (5] which is fitted by a log-normal distribution 
with a mean of log Liso ~ 51.54 (Liso in erg s~*), a dispersion 
o-iogLi.o = 0.795, and x? = 1-24. 

The Swift trigger is very complex and the sensitivity of 
the detector is very diffi cult if not impossible to parameter- 
ize exactly (|Bandl I2OO6I ). However, an effecti ve luminosity 
thres hold appears to be present in the data ()Kistler et al.l 
'2OO8!, their figure 2). I find that the luminosity thresh- 
old can be approximated by a bolometric energy fiux limit 



Llim -- 

then 



1.2 X 10 erg cm s . The luminosity threshold is 



(9) 



where Dl is the luminosity distance to the burst. 

With the above luminosity threshold and an adopted 
GRB rate history, the observed luminosity distribution can 
be fitted by an intrinsic Schechter LF with a power-law index 
— 1.225 and a characteristic luminosity 0.986 x 10^^ erg s~* 
(the dashed line in Fig. (5] for details see Section [6. 2p . 

The data have a moderate excess around Liso = 3.8 x 
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lO'^^erg s~^, at the 2A-a level (relative to the log-normal 
distribution). It probably indicates the existence of a faint 
population of GRBs, which will be discussed in details in 
Section |6l 



4.2 Selection Biases 

Selection effects involved in a GRB sample are hard to 
model quantitatively. There are at least two kinds of se- 
lection effects at work, which have been extensively dis- 
cussed in the lit e rature (iBloom et al.l 20031: JCoward 20071: 
Fiore et a l' '2007: 'Guctta fc Delia Valle"200?: i Kistler et al l 
20081 : iLe fc Do rmer 2007.; Co ward et al. 200^: (1) 



] 

GRB 

detection and localization; and (2) redshift determination 
through spectroscopy /photometry of the GRB afterglow or 
the host galaxy. A full address of these issues is beyond the 
scope of this paper, so I only give a brief description on 
them. 

Since all the GRB s in the sample used in this paper (and 
in lKistler et al.ir2008l ) are detected by Burst Alert Telescope 
(BAT) onboard Swift, biases arising from GRB detection are 
minimized although the sensitivity o f BAT is very difficult 
to parameterize exactly (|Bandll2006l ). An effective luminos - 
ity threshold introduced above and in iKistler et all (|2008l ) 
would be a reasonable approximation of the detection crite- 
ria for the bursts in the sample. 

Due to the transient nature of GRBs, a fast and accu- 
rate location of the burst is crucial for the determination 
of the redshift. The automatic slew ability and the muti- 
wavelength nature of Swift have led that an accurate local- 
ization of the GRB is possible shortly after its detection, 
with an efficiency higher than any previous GRB obser- 
vatory (e.g., BeppoSAX and HETE2). An impressive num- 
ber of ground-based facilities (dedicated robotic telescopes, 
VLT, Gemini, Keck, etc) have been involved in Swift follow- 
up observations. However, optical afterglows have been dis- 
covered for only ~ 50% of Swift GRBs, a fraction only 
sligh tly greater than that of BeppoSAX and HETE2 sam - 
ples (|Berger et a"Lll2005l : lFiore et al.ll2007l : lKann et al]|2008l ). 
The cause for the non-detection of optical afterglows f or well 
localized GRBs remains a puzzle ((Castro- Tirado et al ] |2007l : 
iRol et al.|[2007l : iKann et allboOSD . 

Generally, an apparently bright burst would be easier 
to localize than an apparently faint burst. On average more 
distant bursts are expected to have lower observed gamma- 
ray fluxes, hence the number of GRBs with hight redshift 
might have been somewhat underestimated. But this indi- 
cates that the deviation of the GRB rate from that predicted 
by the SFR could be mor e serious than that has been ob- 
served (|Kistler et al.ll2008l ). 

For the GRBs detected and localized by Swift, most 
of them are foUowed-up spectroscopically by ground-based 
telescopes. The participation of ground-based telescopes in 
follow-up observations has greatly increased the number of 
GRBs with redshifts (about 1/3 of the total Swift GRBs). 
However, this also makes the selection biases in the red- 
shift determination extremely complex, since instrumental 
selection biases of different telescopes have to be involved. 
Different techniques have also been applied in the redshift 
measurement: absorption lines, emission lines and photom- 
etry of the Lyman edge, which causes additional biases. 

In addition to the various selection biases in redshift 



measurements discussed above, there has been a well-known 
'redshift desert' in the redshift interval 1.4 < z < 2.5 
since galaxies in that range have been hard to detect 
spectrosc opically with t raditional ground-b ased telescopes 
(Adelberger et al.l 12004 ISteidel et al]|2004 ). However, the 
'redshift desert' problem does not seem to be serious for 
Swift GRBs, since among the 85 long GRBs with measured 
redshifts by 31 March 2008, 23 bursts are in the redshift 
interval 1.4 < z < 2.5 (Fig. Ilip . This is probably caused 
by the fact that different telescopes have different spectral 
ranges and they complement each other to some degrees. 



5 A MODEL FOR THE GRB RATE HISTORY 



There is evidence that GRBs are beamed (Harrison et al.l 
1 19991 : iKulkarni et"al]|l999l : IStanek et aHliggd '). Hence, when 
discussing the probability distribution function of GRBs, the 
effect of jet beaming must be taken into account. Assuming 
that a GRB radiates its energy into two oppositely directed 
jets, each having a half-opening angle Sjot. The total solid 
angle spanned by the jets is then Ahlj, where the beaming 
factor u = 1 ~ cosSjct < 1. 

For simplicity, I assume that except the comoving rate 
density, the property of GRBs does not evolve with the cos- 
mic redshift. Then, the intrinsic distribution function of z, 
Lisa and y = log (tanSjct) must have a form 



P {z, Liso, y) = f{z)(t> (Liso) ^ {Liso,y) 



where 



jGRB 



(z) rfKc 



(10) 



(11) 



1 + z dz 

(c///o)~'^Sgrb(z) is the comoving rate density of GRBs, 
and Vcoui is the comoving volume. Here c is the speed of 
light. 

In equation (|10p . tp{Liao,y) is normalized with respect 

to y 



ip(Liso,y)dy = 1 



(12) 



In a flat universe (Sim + = 1), the comoving volume 
is calculated by 



com 

J ^com J 1 

dz dz 
where the comoving distance 

^ , dz' 

-L-'com \ ^) 



Ho Jo ^an(l + 2;')'^ + fiA 



(13) 



(14) 



A GRB is detected by an observer only if one of its jets 
points toward the observer. Hence, without consideration 
of other selection effects, the probability for a GRB to be 
detected by an observer is equal to uj. Multiplying equation 
(|10|) by ui then integrating it over y, one gets the observed 
distribution function of z and Lisa without considering other 
selection effects 



P(z,Li.o) =/(z)$(Li,o) , 
where the beaming-convolved LF 

<l>(Liso) = {uj)(t){Liso) , {(^) : 



(15) 



LJXp{Liso,y) dy . (16) 
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threshold for the production of GRBs. Studies on the GRB 
progenitors and the collapsar model predict that long- 
duration GRBs are produced only f o r progenitor stars with 



Z/Zq < 0.1 (Hirschi et all l2005l: lYoon fc Langc"? 



IWooslev fc Heger 2006, ; lYoon et all l2006l : IVink 2003). In 



20051 : 



vestigations on the host galaxies of long GRBs show that 
an upper cut-off for the gas-phase metallicity of GRB hosts 
is likely 12 -f- log(0/H) ~ 8.5, corresponding to Z ^ 0.2- 
0.6.^0 depending on the adopted metallicity scale and solar 
abundance vale (Modjaz et al. 2008; see also Wolf & Podsi- 
adlowski 2007, Savaglio et al. 2008). 

In general, the stellar metallicity is always lower than 
the gas-phase metallicity ([Gallazzi et a l. 2005). Hence the 
two metallicity thresholds from the theoretical study on 
GRB progenitors and the observational investigation on 
GRB hosts appear to be consistent with each other. Since 
the metallicity evolution adopted in equation [T5] and dis- 
cussed in Section [3] refers to the gas-phase metallicity, 
throughout the p aper I assume that e = 0.3 . 

According to lLanger &: NormanI (|2006l ). 



Figure 6. The cosmic volume coordinate Q (defined by eq. [TtJ as 
a function of the cosmic redshift z (solid curve) , and the derivative 
dQ/dz (dashed curve). (Qm = 0.3, Q\ = 0.7.) 



Generally, the luminosity-averaged jet beaming factor 
{uj) is a function of Liso. Therefore, the beaming-convolved 
LF <E> (I/iso) differs from the intrinsic LF tf) (Liso). 

For the purpose of studying the rate density history of 
GRBs and the star formation, it is more convenient to use a 
dimensionless volume coordinate Q than to use the redshift 
z, where Q = Q{z) is defined by l|Kistler et al.ll2008l ') 



1 rfKo, 

1 + z' dz' 



-dz' 



(17) 



which increases monotonically with z (Fig. [6]). For z <C 1, 
one has Q ~ iirz^/S. As 2: — > 00, one has Q constant. 

The coordinate Q is particularly useful in binning the 
data since the definition of Q has taken into account both 
the effect of comoving volume and the effect of cosmic time 
dilation. For example, if the comoving rate density of GRBs 
were a constant, in each equally-sized bin of Q the number 
of GRBs would be a constant. In contrast, if the data are 
binned with an equal size of Az, the number of GRBs would 
change dramatically from bin to bin because of the comoving 
volume and the cosmic dilation factors in equation (|lip . 

By equations (|lip and (|15p . the distribution function of 
Q and Liso is 



dz 

P{Q,L,..) = P[z,L,,^) — 



T.grb{Q)$ {Li, 



(18) 



Following iLanger fc NormanI l|2006l ). I assume that the 
GRB rate is related to the SFR and the gas-phase metallicity 
in the host galaxy by 



Sgrb(z) = A-^{z,e)T.sFR{z) 



(19) 



where Esfr(2) is the comoving rate density of star for- 
mation, ^{z,e) is the fractional mass density belonging to 
metallicities below Z = eZQ at a given redshift z, and A is 
a normalization factor. 

The parameter e is determined by the metallicity 



^{z,e) 



r(Q -f 2) 



(20) 



where a ~ —1.16 is the power- law index in the Schechter 
distribution function of galaxy stellar masses (Panter et al.l 
|2004 ). and /3 « 2 is the slope in the linear bisector fit to 
the galaxy stellar mass-metaU icity relation (|Savaglio et al.l 
I2OO5I ; iLanger fc NormanI l2006h . In equation it is as- 

sumed that the average cosmic metal licity evolve s with red- 
shift lw_-0. 15 dex per unit redshift l|Kewlev fc K obuhuckvl 
l2005l . l2007l . see eq. [7] in Section |3] and the relevant discus- 
sion) . 

For the SFR. lKistler et all ll2008t) adopted t he piecewise 
power- law derived by Hopkins fc BeacomI (|2006l ). Here I use 
the the piecewise power-law derived from the updated data, 
given by equations ^ and ([2]). The two formulas differ sig- 
nificantly only at 2: > 4, where measurements of the SFR 
are highly uncertain. 

For the beaming-convolved LF of GRBs, I assume that 
it has a Schechter-function form 



$ (Liso) 



exp (-Liso/L*) 



(21) 



where ql and are constant parameters to be determined 
by observational data. Although people often adopt a bro- 
ken power-law or a double power-law f or the GRB LF (e.g., 
iGuetta et aP l2005l : iLiang et al.l l2007l l. I find that a sim- 
pler Schechter function is enoughQ Indeed, it appears that 
a power-law GRB LF at high end tends to over-produce 
the number of bright GRBs (B. Zhang, private communica- 
tions). 



iNataraian et al.l (jSooi) and ISalvaterra fc Chincarinil (|2007|) 

also assumed a Schechter function for the GRB luminosity in 
their work. 
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6 FITTING THE REDSHIFT DISTRIBUTION 
OF SWIFT GRBS 

6.1 Bright GRBs 

To avoid the complication of a detailed treatment of the 
Swtft det ector's threshold a nd the assumption about the 
GRB LF. lKistler et alj (|2008f ) adopted a model-independent 
approach by selecting only GRBs with Liso > 10^^ erg s~^ 
and 2 < 4. The cut in luminosity and redshift minimizes the 
selection effect in the GR B data. With 36 Swift GRBs that 
satisfy the above criteria, iKistler et al.l ((2008) showed that 
the rate of GRBs increases with the redshift much faster 
than the SFR. A Kolmogorov-Smirnov test revealed that 
the SFR alone is inconsistent with the GRB rate history at 
the 95% level. 

For GRBs with Liso > Lum, the observed number of 
GRBs in an observer's time interval Aiobs and with Q in 
the interval Q-{Q + dQ) is N{Q)AtohB dQ, where 

N{Q) = / P{Q,Li,a)dUso 



(22) 



= Sgrb(O) / $(iiso)dLiso . 
Submitting equation (|21|l into equation (|22|) . one gets 



iV(Q) = Sgrb 



(Q)r( 



1 + "L, 



(23) 



There are 45 GRBs with Liso > lO'^^erg s~^ in the Swift 
sample, versus the 44 GRBs in iKistler et all (120081 ). The 
difference in the two numbers is caused by GRB 050318, 
whose luminosity is very close to lO^^erg s~^. IKistler et al.l 
l|2008h used the isotropic- e quival ent energy of GRB 050318 
published in iButler et al.l ([2003) to calculate the luminos- 
ity and obtained a value that is slightly below the luminos- 
ity cut. As mentioned in Sectio n | 4. II I take the isotropic - 
equivalent energy of GRBs from lKocevski fc Butlej (|200'i1 ). 
I get a luminosity for GRB 050318 that is slightly larger 
than lO^^erg s~^. To avoid this ambiguity arising from the 
luminosity uncertainty, hereafter I assume a luminosity cut 
Liim = 0.8 X lO^^erg s~^ Then, the total number of GRBs 
with Liso > Liini is unambiguously 45. 

The distribution of Q for the Swift 45 GRBs is plotted 
in Fig. [3 The additional GRB 050318 falls in the third bin, 
resulting that the nu mber of GRB s in th e third bin is larger 
by one than that in IKistler et al.l (|2008l ) (compare to their 
fig. 3). 

Since Lum is a constant, by equation (|22p the number 
of GRBs in each equally-sized bin of Q is proportional to 
Sgrb, independent of the GRB LF. For the selected bright 
GRBs, the luminosity threshold a rising from the det ector 
flux limit has no effect when z<A l|Kistler et al.ll2008l '). Us- 
ing the Sgrb given by equation (|19|) with Ssfr given by 
equations ([TJ and ^ (and q = — 1.16, /3 = 2) to fit the first 
six data points in Fig. [7] which have 2 < 4 by varying only 
the normalization, I get a surprisingly good fit as shown by 
the dashed curve in the figure, with Xi ~ 0.14. This fact in- 
dicates that the GRB rate density assumed in equation (I19|) . 
which takes into account the evolution of the cosmic metal- 
licity, reasonably represents the true GRB rate. Hence, I will 
adopt this GRB rate density in all the following analysis. 



15 



10 



Figure 7. Distribution of Q for 45 Swift GRBs with Liso > 
0.8 X lO^'^erg s~^ (the solid histogram, with the number of GRBs 
is each bin indicated by a dark point with Poisson error bars). The 
dashed curve is the best fit of the Sgrb (eg. 1191 with EsFR given 
by eqs. [T]and[2ll to the first 6 data points {z < 4) by varying the 
normaUzation, which has = 0.14. The deviation of the model 
from the data at 2 > 4 is presumably caused by the flux limit of 
the detector which results a decrease in the number of detected 
GRBs. 



6.2 All GRBs 

Now 1 consider the effect of the detector flux limit (or, equiv- 
alently, the luminosity threshold) on the distribution of lu- 
minosity and redshift for the 64 GRBs in the sample. As 
mentioned in Section [4. II the observed luminosity threshold 
can be modeled by an energy flux limit Fnm ~ 1.2 x 10~*erg 
cm^'^ s^^. Then, Lum is a function of z (eq.[9]). The observed 
distribution of Liso is then given by 



^ (Liso) = $ (Liso) Atobs 



f{z)dz , 



(24) 



where 2max = Zmax (Liso) IS the maximum redshift up to 
which a GRB of luminosity Liso can be detected, solved from 
equation Liim(2) = Liso. 

Submitting equation (|21|l into equation ([24} then fitting 
the data in Fig. [Sjwith $, I get ol = -1.225, L* = 0.986 x 
lO'^^erg s-\ A' = AtobsA = 9.386, and x? = 1-18 (the 
dashed curve in Fig. [5]). The data excess around Liso ~ 3.8 x 
10*^ erg s~^ is at the 2.3-0- level relative to the modellf] 

The N{Q) is still given by equation (|23|l . but now Liim 
is a function of Q. Since —2 < a < —1, the incomplete 
gamma function in equation (|23p can be evaluated by the 



^ I tried also to fit the data without the outlier point. But the 
main results in the paper are not affected drastically. 
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Figure 8. Distribution of Q for all the 64 Swi/t GRBs (the solid 
histogram, with the number of GRBs in each bin indicated by a 
dark point with Poisson error bars). The dashed curve shows the 
distribution given by equation 1231 1 with the parameters obtained 
by fitting the luminosity distribution (Fig. [5] the dashed curve). 
The first data point has an offset about 2.7-a from the dashed 
curve. The overall = 1-01. If the first data point is not included 
in the calculation of chi-squares, then X? = 0.24. 



recurrence formula 

r(a,a;) = i fr(l + a,a;)-a;'"e-"l . (25) 

In Fig. |8] I show the distribution of Q for all the 64 
Swift GRBs in the sample. The dashed curve is the N(Q) 
calculated by equation (|23p with the normalization and the 
LF parameters determined above, and Lum calculated by 
equation ((9)|. Globally, the modeled N{Q) fits the observa- 
tional data very well (without adjustment of parameters), 
with Xi ~ 1-01. However, there is an obvious excess in the 
number of GRBs in the bin of < Q < 1, at the 2.7-a level. 
If the first bin is excluded in calculating the chi-squares, I 
get Xr = 0.24. 

The rate of bright GRBs in equation (|23|) with the nor- 
malization obtained above and a constant Lum = 0.8 x 
lO^^erg s~^ fits the first six data points {z < 4) in Fig. [7] 
with x'^/d.o.f — 0.682/6, very close to the best fit ob- 
tained by varying the normalization in Fig. [7] which has 
xVd-o.f = 0.681/5. 



6.3 The Cumulative Distribution of Q 

The cumulative distribution of Q for 32 Smft GRBs with 
Liso > 0.8 X lO^^erg s"^ and z < 3.5 is shown in Fig. [90 



° I set the upper bound of the redshift to 3.5 instead of 4 to 
reduce the effect of the luminosity threshold. 



z 



0.5 1 1.5 2 2.5 3 3.5 




12 3 4 



Q 

Figure 9. Cumulative distribution of Q{z) for 32 Swift GRBs 
with z < 3.5 and Liso > ^min = 0-8 x lO^^erg s~^ (the stepwise 
curve). The cut in redshift and luminosity is chosen so that the 
data are not affected strongly by the luminosity threshold. The 
smooth solid curve is calculated with the GRB rate Sgrb in 
equation I I19I I. which is independent of the GRB LF. The dashed 
curve is calculated with the SFR given by equations ((TJ and Q 
alone. 

This subsample of GRBs is not subject to the luminosity 
threshold. The integral distribution of Q, defined by 

N{<z)= / N{Q)dQ , (26) 
Jo 

is shown with a smooth solid curve, which agrees with the 
observed data very well. The integral distribution of Q given 
by the SFR alone (the dashed curve) does not agree with the 
observation. 

In Fig. 1101 I show the cumulative distribution of Q for 
all the 64 GRBs in the sample. To show the excess of GRBs 
at low redshift, I use the cumulative distribution defined by 

N{>z)= / N{Q)dQ, (27) 

Jq(z) 

where Zmax = 7. The N{> z) given by the model is shown 
with a smooth solid curve, which fits the observed distribu- 
tion beyond z — 0.7 very well. An excess in the number of 
GRBs at redshift < 0.7 is clearly seen. 

Given the small number of GRBs in the sample, the ob- 
served excess in the number of GRBs at low redshift (and low 
luminosity, Fig. [5} could simply arise from statistical fluctu- 
ations. However, it is also possible that the excess reflects a 
real deviation of the redshift (and luminosity) distribution 
from the model, because of the following facts: 
(A) The detection of highly sub-luminous and sub-energetic 
nearby GRBs 980425, 031203 and 060218 has led people 
to propose that there exists a unique population of faint 
GRBs, whose rate is much higher than normal cosmological 
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Figure 10. Cumulative distribution of Q(z) for all the 64 GRBs 
in the sample (the stepwise curve), with N{> z) defined by equa- 
tion ((27J ( ^max — 7) . The smooth solid curve is the result given by 
the model. The vertical dotted line denotes z = 0.7 {Q = 0.564). 
An excess in the number of GRBs at 2 < 0.7 is clear (c.f. Fig.[8ll. 



GRB s llCobb et al.ll2006al: iPian et ai]|2006l: ISoderberg et all 
20061: IChapman et all l2007l : iGuetta fc Delia Valid l2007l : 



Liang et al 



20071 ) 



(B) Some nearby long-duration GRBs are found not to have 
accompanied supernovae and hence are probably not re- 
lated to the death of massive stars, challenging the tradi- 
tional scheme for the classification of GRBs by their dura - 
tions l|Gehrels et al.ll2006l : [Zhan3l2006l : [zhang et al.ll2007al ). 
These non-supernova GRBs include 060614 at z = 0.125, 
060505 at z = 0.0 89 (Cob b et all l2006bl: [Pella Valle et all 
l2006l : lFvnbo et al.li2006b: ,Gehrels et al.ll2006l ). and 051109B 
at z = 0.08 (|Perlev et al lfiooi T 

(C) The SFH of the Local Group and other nearby galax- 
ies indicates an excess of the local SFR density rela- 
tive to the cosmic SFH in the recent epoch of z < 0.5 
l|Drozdovskv et ai]|2008l ). This might also be a cause for the 
excess in the number of GRBs at z < 0.7. 



7 CONCLUSIONS 

I have presented an updated cosmic SFH up to redshift z — 
7.4. The updated sample of SFR data are obtained by adding 
the new UV and IR measu reme nts on the SFR de nsity of 
iBouwens etHI (|2007l . |2008| ) an dlReddv et al.l (|2008[| to th e 
'good' data sample compiled bv lHopkins fc Beacom (l2006h. 
The two joined UV-I-IR measurements of iReddv et al.l 
l|2008t ) at z ^ 2.3 and z ~ 3.05 agree well with previous 
measurements of SFRs in the redshift interval 1 ^ z < 4, 
and are consistent with a flat evolution SFH i n this inter- 
val. The UV measurements of iBouwens etHI (|2007l , |2008| ) 
at z ~ 3.8, 5, 5.9, 7.4 and an estimate at 2 ~ 10 significantly 



expand the SFR data sample at 2: > 4, and are broadly con- 
sistent with the previous results in 3.5 ^ z < 6. The updated 
sample provides a consistent picture for the cosmic SFH up 
to 2 ~ 7.4 (Fig. [Sjl , although the dust correction at 2: > 4 
is highly uncertain w hich results large uncert ainties in the 
SFH at high redshift l|Drozdovskv et al.ll2008l ). 

The updated SFH still under-produces the GRB rate 
density at high redshift when compared to the Swift GRB 
redshift distribution (Fig. [1]), confirming the previous claim 
toaigne et al.ll2007l : iLe fc Dermeill2007l : lKis"tler et al.ll2008l ) . 
The discrepancy is investigated under the assumption that 
long-duration GRBs trace both the star formation and the 
metallicity evolution, as motivated by the observations that 
long GRBs occurred in star-forming galaxies with low metal- 
licities and the theoretical study on the coUapsar model 
which shows that GRBs can be produced only from low- 
met allicity massive stars. 

Since the cosmic metallicity decreas es with redshift 
llPettini et al. 1999; Prochaska ct al. 2003; Rao et al.|[2003l ; 
Kobulnic kv fc Kewl ev 2004; IKewlev fc Kobulnicky 200- 
2007; Kulkarni et al. .2005. ; iSavaglio et al., .2005,; ,Savagli 
2006), it is natural to expect that the ratio of the GRB 
rate to the SFR must increase with redshift if the sce- 
nario that long GRBs are pr oduced by the death of mas- 
sive stars with low metallicity l|Mac Fadvcn fc Woosley 1999|; 
MacFadven et al.ll200ll: i Hirschi et al..,2005 : .Yoon fc Lange3 



20051 : IWooslev fc Bloomll2006l ) is correct. Adopting a simple 



model for relating the GRB rate densit y to the SFR den- 
sity a nd the cosmic metallicity evolution (|Langer fc Norm"anl 
I2OO6I ) and assuming a flux limit for the Swift detector, I have 
shown that the redshift distribution of the 64 Swift GRBs 
with measured redshifts and calculated luminosities can be 
successfully fitted by the updated SFH with a threshold in 
the metallicity for GRB production (Figs. [T WTU)) . 

iKistler et al.l l|2008l ) have considered several possibilities 
for the cause of the discrepancy between the the Swift GRB 
rate and the SFH. They have shown that the Kolmogorov- 
Smirnov test disfavors an interpretation as a statistical 
anomaly. Selection effects are also not likely to cause an 
increased efficiency in detecting hight- redshift GRBs. Al- 
though Kistler et al. have argued that alternative causes are 
possible (e.g., evolution in the fraction of binary systems, an 
evolving IMF of stars, etc), the results in this paper indicate 
that the cosmic metallicity evolution may be the simplest in- 
terpretation. 

However, the results show an excess in the number of 
GRBs with low luminosity (Liso ~ 3.8 x 10*® erg s~^; Fig.[5j 
and at low redshifts {z < 0.7; Figs. [S] and [TDJ. The existence 
of an excess is confirmed by the up-to-date Swift GRBs with 
measured redshifts, detected by 31 March 2008 (Fig.ITTJ- Al- 
though it might simply be caused by statistical fluctuations, 
the observed excess could also be consistent with the specu- 
lation that there is a unique population of intrinsically faint 
and nearby GRBs (Section [S}. 

The GRB s ample used in this paper (and in 
iKistler et al.ll2008l ) is almost definitely incomplete and non- 
uniform, because of the complex selection biases in the red- 
shift measurement discussed in Section 14.21 In fact, all the 
current works on the redshift distribution of GRBs have such 
a problem. However, IKistler et al.l (|2008l ) have argued that 
none of the selection biases appears to be able to increase 
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Figure 11. Distribution of Q for all the 85 Swift GRBs with red- 
shifts, detected by 31 March 2008 between and including GRBs 
041220 and 080330 (the solid histogram). The dashed curve shows 
the distribution given by equation I I23I I with the parameters ob- 
tained by fitting the luminosity distribution (Fig. |5] the dashed 
curve), renormalized by the total number of GRBs. The first data 
point has an offset about 3.0-0" from the dashed curve. The overall 
Xr = 1.21. If the first point is not included in the calculation of 



chi-squares, then xf 



■ 0.28. 



the overall observability with redshift and account for the 
enhancement in the GRB rate relative to the SFR. 

The number of GRBs in the sample is small, virhich also 
prohibits one from obtaining a strict constraint on the pa- 
rameters in the GRB LF and the cosmic metallicity evolu- 
tion. 

Despite the above problems, the results of this work 
suggest that the rise of the observed Smft GRB rate rel- 
ative to the SFR is compatible with an interpretation by 
the evolution of the cosmic metallicity. Once the problems 
are solved or significantly alleviated in future, a significantly 
improved and enlarged sample of GRBs with measured spec- 
tra and redshifts will become available. Then, by comparing 
the observed GRB rate history to the SFH determined with 
other approaches, it will be possible to probe the cosmic 
metallicity evolution with GRBs. 
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